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We study the ground state of the one-dimensional half-filled Hubbard model with on-site (nearest- 
neighbor) repulsive interaction U {V) and nearest-neighbor hopping t. In order to obtain an accurate 
phase diagram, we consider various physical quantities such as the charge gap, spin gap, Luttinger- 
liquid exponents, and bond-order-wave (BOW) order parameter using the density-matrix renormal- 
ization group technique. We confirm that the BOW phase appears in a substantial region between 
the charge- density- wave (CDW) and spin- density- wave phases. Each phase boundary is determined 
by multiple means and it allows us to do a cross-check to demonstrate the validity of our estimations. 
Thus, our results agree quantitatively with the renormalization group results in the weak-coupling 
regime {U ^ 2t), with the perturbation results in the strong-coupling regime {U ^ 6t), and with the 
quantum Monte Carlo results in the intermediate-coupling regime. We also find that the BOW-CDW 
transition changes from continuous to first order at the tricritical point {Ut^Vt) ~ (5.89t, 3.10t) and 
the BOW phase vanishes at the critical end point {Uc^Vc) ~ (9.25t, 4.76t). 

PACS numbers: 71.10.Fd, 71.10.Pm, 71.10.Pm, 71.30.+h 



For several decades quasi-one-dimensional (ID) ma- 
terials, e.g., organic conductors [l|, conjugated poly- 
mers [2|, and carbon nanotubes [3|, have been a main 
subject of research in the field of condensed matter 
physics. A minimal electronic model which can describe 
their basic properties is the ID extended Hubbard model 
(EHM) [3] . The Hamiltonian is given by 

+ U^ni^riii ^V^niaUi^ia^, (1) 



where cj^ (Qct) is creation (annihilation) operator of an 
electron with spin a at site and n^cr = 4a number 
operator, t is nearest-neighbor hopping term and U (V) 
is on-site (nearest-neighbor) Coulomb interaction. De- 
spite the geometric simplicity, this model at half filling 
is believed to exhibit a variety of phases due to strong 
quantum fluctuations. 

Within the g-ology scheme the system has merely 
two insulating phases when the interaction strengths 
are positive: for U < 2V the ground state is 2/cf- 
charge- density- wave (CDW), where both the charge and 
spin excitations are gapped; for U > 2V a Mott in- 
sulator with 2/cF-spin-density-wave (SDW), where the 
spin excitation has no gap. However, based on non- 
perturbative numerical results, Nakamura argued that 
there is also a bond-order- wave (BOW) phase, where the 
ground state has a long-range staggered bond order, be- 
tween the CDW and SDW phases So far much ef- 
fort has been devoted to fix the ground-state phase dia- 
gram both analytically 0, H, 0, [HI, [lH and numeri- 
cally H, [3, [H, Nevertheless, surprisingly their 
results are in few (quantitative) agreements with each 
other. The aim of this paper is to produce a highly ac- 



curate phase diagram of the ID half-filled EHM and to 
resolve the apparent contradictions. 

We employ the density-matrix renormalization group 
(DMRG) method, which is one of the most powerful 
numerical techniques for studying ID many-body sys- 
tems [l8|. With open-end boundary conditions, ground- 
state and low-lying excited-states energies as well as ex- 
pectation values of physical quantities can be obtained 
quite accurately for very large finite-size systems (up 
to sites L - O(IOOO)). In DMRG procedure we keep 
m = 1200 to 3000 density-matrix eigenstates, which are 
much larger than those in the previous DMRG stud- 
ies [li, and all the calculated quantities are ex- 
trapolated to the m ^ oo limit. In this way, the maxi- 
mum truncation error, i.e., the discarded weight, is less 
than 1 X 10~^^, while the maximum error in the ground- 
state energy is AE/t ~ 10~^ — lO"*". We strongly argue 
that such large m values and the m-extrapolation are 
essential for required accuracy of the measurements. 

In order to determine the phase diagram including two 
phase boundaries, we calculate several physical quanti- 
ties. Each boundary is determined by multiple means 
from the quantities and it allows us to do a cross-check on 
the estimates. First, to obtain the BOW-CDW boundary 
we calculate the charge gap 

Ac= lim [£;(A/' + 2,0) + £;(7V-2,0)-2£;(7V,0)]/2, (2) 

L^oo 

where E{Nq^Sz) is the ground-state energy for a given 
number of electrons Nq and z-component of total spin 
Sz' We take N = L for half- filled case. In the atomic 
limit t = 0, the phase boundary becomes a line U = 2V 
with Ac = U{= 2V). If finite t is introduced, the system 
can gain some kinetic energy of the order of t near the 
BOW-CDW instability due to the competition between 
the on-site and nearest-neighbor Coulomb interactions. 
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Thus, the charge gap is minimized at the BOW-CDW 
boundary. Next, to evaluate the SDW-BOW boundary 
we calculate the spin gap 



lim [E{N,1) - E{N,0)]. 



(3) 



If y <C the system is a Mott insulator with 2/cf- 

SDW. The electrons are uniformly distributed over the 
system, so that there is no spin gap. As V increases, 
the charge fluctuations are enhanced, and then a transi- 
tion from the SDW phase to the BOW phase occurs. In 
the BOW phase, the electrons polarize alternatively and 
spin-singlet bound states are formed on dimers. Con- 
sequently, we can make an estimate of the SDW-BOW 
boundary as a point where the spin gap begins to de- 
velop. However, for some parameters the spin gap is too 
small to figure out if it remains finite, i.e., Ag ^ 10~^t. 
Therefore, for verifying the presence of the spin gap we 
consider the spin-spin correlation function 



Siq) 



Mik-l) 



iisls^) - (4) (s!)) (4) 
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with q = 27r/L and = — n^|. According to the 
Luttinger liquid theory [3], the long-range behavior of 
this function is governed by the LL spin exponents K^j 
(= limq^o 7TS{q)/q). We find K^j = in the spin-gapped 
phase and 1 everywhere else in the thermody- 

namic limit [2^. This criterion enables us to estimate 
the SDW-BOW critical point precisely. Although we can 
obtain all the phase boundaries with the quantities men- 
tioned above, the BOW oder parameter is also studied 
for making extra sure. The order parameter simply gives 
the boundaries between the BOW phase and the other 
phases. The BOW operator is given as 



0. 



(5) 



and we define the BOW order parameter {B) as an ampli- 
tude of the BOW oscillation in the center of the system, 



i.e., (B) = limL^c 



(B 



L/2 



5l/2+i)|. For {B) ^ 0, a 



long-range order of the BOW state appears. 

A careful extrapolation of these quantities is neces- 
sary to extract correct value in the thermodynamic limit 
L ^ oo. We thus study various lengths of chains with 
L = 32 to 512 and perform finite-size-scaling analysis 
based on the L-dependence of the quantities. Figure [T] 
shows the finite-size-scaling analyses for (a) the charge 
gap, (b) spin gap, (c) spin-spin correlation function, and 
(d) BOW parameter near the phase transitions at = 4. 
The charge (spin) gap is systematically extrapolated by 
performing a least-squares fit to the forth-order polyno- 
mial in refiecting the holon (spinon) band structure 
around the band edge. Then, an estimation of the LL 
spin exponent in the thermodynamic limit is not so sim- 
ple for finite-size calculations. In the spin-gapless phase. 
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FIG. 1: Finite-size-scaling analyses for (a) the charge gap, 
(b) spin gap, (c) spin-spin correlation function, and (d) BOW 
order parameter near the phase boundaries at U — 4t. 



one cannot expect easily find Kcr 1 exactly due to log- 
arithmic corrections. However, the logarithmic correc- 
tions are known to vanish at which the spin gap opens, 
in analogy with the dimerization transition in the J1 — J2 
model [21]. In the spin-gapped phase, there is a similar 
difficulty as follows; if the spin gap is small, the conver- 
gence of i^cr to will obviously occur only for very large 
systems. As a result, we will estimate the critical point 
where the spin gap opens by 7TS{q)/q crossing 1 at g ^ 0. 
This method was primarily used in Ref. [14]. Let us now 
turn to the BOW order parameter. Since the order pa- 
rameter in the thermodynamic limit is very small com- 
pared to the finite-size results, a well-controlled finite-size 
extrapolation is mandatory. In our calculations, the most 
problematic finite-size effects are the Friedel oscillation 
due to the open edges. Assuming that the amplitude of 
the Friedel oscillation in the center of a finite chain scales 
as L~^p [22[, the BOW order parameter would be well- 
extrapolated as a function of For example, we 
may expect Kp ^ 0.5 in the vicinity of the SDW phase, 
so that (B) is scaled better by 1/vT than by 1/L near 
the SDW phase. 

Figure [2] shows the extrapolated results of (a) the 
charge gap, (b) spin gap, (c) spin correlation function, 
and (d) BOW order parameter around the phase transi- 
tions {U ~ 2V) as a function of V/t for U = 4t. Let us 
look at the charge gap to estimate the BOW-CDW phase 
boundary. The charge gap decreases with approaching 
to a point V ~ 2.164t and vanishes smoothly at the 
point. In other words, both the BOW and CDW insulat- 
ing gaps start to develop gradually at the point. It means 
that a continuous transition between the BOW and CDW 
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FIG. 2: Extrapolated results of (a) the charge gap, (b) spin 
gap, (c) spin-spin correlation function, and (d) BOW order 
parameter near the phase transition for [/ = 4t. The dashed 
and dotted lines denote the SDW-BOW and BOW-CDW crit- 
ical points, respectively. Insets: same quantities plotted with 
another scale. 



phases occurs at the critical point V ~ 2.164t. Note that 
the BOW insulating gap is of the nature of the Mott type. 
We now turn to the SDW-BOW phase boundary. We find 
that the spin gap is finite for F ^ [7/2 and decreases with 
decreasing V. The critical point appears to lie around 

V = 1.9t from the disappearance point of the spin gap. 
The crossing point with Trlinig^o S{q)/q = 1 gives more 
precise estimation of the critical point V ^ 1.877t. Cor- 
respondent ly, the BOW order parameter has finite val- 
ues only in the region 1.877t < V < 2.164t. With in- 
creasing y, {B) rises exponentially from the SDW-BOW 
critical point, reaches the maximum value ~ 0.18 around 

V = 2.14t, and goes down to zero at the BOW-CDW crit- 
ical point. Note that both values of the critical points are 
in good agreement with those of the previous quantum 
Monte Carlo (QMC) study [H. 

Figure [3] shows the same quantities as in Fig. [2] but 
for U = 8t. Near the SDW-BOW phase boundary 

V ^ 4.039t, the behavior of all the quantities is quali- 
tatively similar to those in the case of C/ = 4t. On the 
other hand, the physical properties seem to be discon- 
tinuous at the BOW-CDW phase boundary V ^ 4.142t, 
which indicates that the transition is of first order. At 
the boundary, the charge gap remains finite and the slope 
of Ac with respect to V is discontinuous. However, the 
value of Ac must be continuous since a competition be- 
tween two kinds of charge configuration, i.e., CDW and 
uniform, leads to the BOW-CDW transition. Associated 
with this charge redistribution, the spin gap jumps by 
two orders of magnitude. In the CDW phase, it comes 
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FIG. 3: The same quantities as in Fig.[2]but for U = St. Solid 
Une in the inset of (b) denotes the spin gap in the V/U ^ 00 
limit, i.e.. As = 3V - U. 



rapidly close to a line Ag = 3^ — U which becomes exact 
in the V/U ^ 00 limit. Also, the BOW order parame- 
ter develops with approaching the BOW-CDW boundary 
and disappears at the transition point. 

Whereas the BOW-CDW transition is continuous for 
U = 4t, it is of first order for U = St. Hence, a tricrit- 
ical point {Ut, Vt), at which the transition changes from 
continuous to first order, must exist on the BOW-CDW 
boundary, as suggested in Refs. [ii,!!!. To evaluate the 
tricritical point, we examine the LL charge exponent Ko 
via the derivative of charge structure factor at g = [23| 



lim - ( 

ki 



-^\{nkni) - {uk) {ni)) . (6) 



Note that Kp is finite only in the continuous Gaussian 
critical point [ll| for small U and zero everywhere 
else. It was shown that the LL exponents can be obtained 
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FIG. 4: Extrapolated results of the LL charge exponent (a) 
and the BOW order parameter (b) on the BOW-CDW bound- 
ary line. Inset: expanded view around the tricritical point 
Ut = 5.89t. 
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quite accurately with DMRG method In Fig. [l^a), 
we plot DMRG results of Kp as a function of U/t on 
the BOW-CDW boundary line. As U/t increases, Kp de- 
creases from 1, reaches to 1/4 at (/7t, Vt) = (5.89t, S.lOt), 
and drops discontinuously to 0; namely, a metal-insulator 
transition occurs at [/ = Ut. Moreover, the Kp curve is 
well-fitted by a function Kp-l/4 = Oml^/{Ut - U)/t 
near the tricritical point [see inset of Fig.[Ha)]. It implies 
that the transition is of the Kosterlitz-Thouless type. 
Let us now consider a point at which the BOW phase 
shrinks to 0, which is called a "critical end point". The 
BOW state is still stable around the tricritical point and 
therefore the critical end point (t/cK) would exist for 
> Ut. For a fixed U (> Ut), the BOW order pa- 
rameter has a maximum around the BOW-CDW bound- 
ary. To find the critical end point, we plot (B) on the 
BOW-CDW boundary as a function of U/t in Fig.H^b). 
{B) decreases with increasing U/t and reaches to at 
{Uc,Vc) = (9.25t,4.76t). For U > Uc, the transition is 
always first-order SDW-CDW one. 




U/t 

FIG. 5: DMRG phase diagram of the ID half-filled EHM. 
The BOW phase exists between the SDW and CDW phases. 

In Fig. [5] we sum up our results as the ground-state 
phase diagram. One can see good agreement with the 
weak-coupling renormalization group (RG) results ^] as 
well as the strong-coupling perturbation results . The 
BOW phase has a maximum width at /7 ~ 4t, which 
is concerned with the fact that the effective nearest- 
neighbor exchange interaction is the largest at the in- 
termediate couplings of U in the half-filled Hubbard 
model It is so because the large exchange inter- 

action promotes the formation of spin-singlet pair if the 
charge fluctuation is introduced by V. Accordingly, we 
confirm that the magnitude of the spin gap is maximized 
around U ~ 4t in the BOW phase. 

In summary, we study the ground-state phase diagram 
of the ID half-filled EHM using DMRG method. We cal- 



culate several quantities with considerable accuracy to 
determine the SDW-BOW and BOW-CDW boundaries. 
As for the phase boundaries, our data agrees quantita- 
tively with the RG results in the weak-coupling regime 
{U ^ 2t), with the perturbation results in the strong- 
coupling regime {U ^ 6t), and with the QMC results 
in the intermediate-coupling regime. We also find that 
the BOW-CDW transition changes from continuous to 
first order at the tricritical point {Ut, Vt) = (5.89t, 3.10t) 
and it locates far from the critical end point {Uc,Vc) = 
(9.25t, 4.76t). Since the previous DMRG results could be 
insufficient in accuracy, our results are not in agreement 
with them. We thus believe that our DMRG results bring 
a sound conclusion and put an end to the controversy on 
the phase diagram of the ID half- filled EHM. 
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